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Abstract. Parker instability leads to the formation of tangential discontinuities in a magnetic field and subsequent 
magnetic reconnection due to a numerical and/or an explicit resistivity. In this paper we investigate the role of 
the uniform, localized and numerical resistivity on the diffusion of magnetic field lines during the growth phase 
of Parker instability modes. We propose a new method to quantify the diffusion of magnetic field lines which is 
attributed to the presence of resistivity in ideal and non-ideal MHD codes. The method relies (1) on integration of 
magnetic lines in between periodic boundaries, (2) on measurements of the dispersion of magnetic field lines with 
the left and the right periodic boundaries and (3) on a statistical analysis of shifts of a large set of magnetic lines. 
The proposed method makes it possible to detect topological evolution of magnetic field. We perform a series of 
resistive MHD simulations of the Parker instability in uniformly rotating galactic disks. We follow the topological 
evolution of the magnetic field evolving due to the Parker instability and relate it to the ratio of total to uniform 
magnetic field in galactic disks. We find that after the onset of the Parker instability, the magnetic field becomes 
first tangled and later on it evolves toward a uniform state due to the presence of resistivity. A similar efi'ect of a 
varying contribution of a turbulent magnetic field is observed in arms and inter-arm regions of galaxies. 

Key words. ISM: magnetic diffusion - Galaxies: structure - Magnetohydrodynamics (MHD) 



1. Introduction 

The issue of stability of magnetized astrophysical disks 
is of primary importance for their dynamical evolution. 
Parker H1966L 11967(1 discovered that galactic disks con- 
taining a large scale horizontal magnetic field and ver- 
tical disk gravity are unstable against the buoyancy of 
magnetic field and cosmic rays. When the contribution of 
the weightless magnetic field and cosmic rays to the total 
pressure is significant, small perturbations of the system 
grow in amplitude, form wavelike perturbations of mag- 
netic field lines and promote sliding of the heavy gas to 
magnetic valleys. The typical horizontal wavelength of un- 
stable modes is of the order of a few hundreds parsecs up 
to a few kiloparsecs. This kind of instability transports 
the magnetic field out of the disks. On the other hand 
Parker instability is supposed to contribute to the mag- 
netic field amplification through the process of fast galac- 
tic dynamo introduced originally by Parker H1992p and 
developed subsequently in a series of papers by Hanasz 
& Lesch [TUM Onnn|) and Hanasz, Otmianowska- 

Mazur & Lesch H2002L hereafter referred to as Paper I), 
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and by Brandenburg and Schmitt H1998II and Moss et al. 

Another relevant MHD instability is the magneto- 
rotational instability, referred to as Balbus-Hawley insta- 
bility H1991inT)92|l . which relies on the angular momentum 
exchange, by means of perturbed magnetic field lines, be- 
tween different orbits of differentially rotating magnetized 
gaseous disks. The development of this instability initiates 
and sustains turbulence and leads to the amplification of 
the small scale magnetic field up to a certain level of sat- 
uration. The common feature of both instabilities is the 
development of non-uniformities and contact discontinu- 
ities on the initially smooth magnetic fields. 

According to Ampere's law, magnetic field discontinu- 
ities correspond to current sheets. Currents are dissipated 
according to Ohm's law and magnetic energy is converted 
into the thermal energy. As a consequence magnetic field 
is no longer frozen in the gas flows and the magnetic 
field topology can evolve due to a nonidealness of the 
medium. Resistive processes may influence the develop- 
ment of MHD instabilities in astrophysical conditions. The 
issue has been undertaken in both cases of Balbus-Hawley 
(Fleming et al. I2000|l and Parker instabilities (Paper I). 

A significant part of studies of MHD instabilities in as- 
trophysics have been done with the aid of numerical sim- 
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ulations. The intrinsic property of numerical codes solv- 
ing MHD equations is the presence of a code-dependent 
numerical diffusivity. Therefore, even if the ideal set of 
equations is actually being solved, the numerical diffusion 
of the magnetic field changes the magnetic field topology. 
In view of the above facts, the main goal of the present 
project is to analyze the influence of the different kinds 
of resistivity on the magnetic field topology in numerical 
simulations of the Parker instability. 

We apply three types of a resistivity: the uniform one, 
the localized resistivity operating only above some criti- 
cal value of the current density, and the numerical one. 
Our model study is based on the general setup similar to 
that applied in Paper I, i.e. we consider a simple equilib- 
rium state in a uniform gravity, assuming that magnetic 
pressure is initially proportional to the gas pressure. The 
initial magnetic field is vertically stratified, and is parallel 
to the azimuthal direction. We superpose a perturbation 
of velocity and simulate the Parker instability in a local 
Cartesian box with coordinates {x, y, z) parallel locally to 
the directions of r, (f> and z coordinates. We take into ac- 
count the disk rotation introducing the Coriolis force. We 
apply the resistivity model following Paper I (and refer- 
ences therein). 

In order to perform a quantitative estimation of the ef- 
fect of resistivity in a Parker unstable system we introduce 
a new method of measuring the topological variations of 
magnetic field lines. The method described in Section O 
allows us to examine the effects of different prescriptions 
of the resistivity and to compare them to the effect of the 
numerical resistivity. The applied concept is similar to the 
description of magnetic diffusion with the aid of Liapunov 
exponents (Barge et al. . 119841 Barghouty and Jokipii 1996', 
Zimbardo ilQOSi) . Our method is different from the one 
applied by Fleming et al. I|2()()()|l . where the resistivity is 
measured through the determination of the saturation of 
the growth of magnetic energy as the result of the Balbus- 
Hawley instability. In the case of Parker instability such a 
method does not provide any information about the resis- 
tivity since without the action of differential rotation we 
do not observe any growth of the magnetic energy. 

We are primarily interested in the role of the resistiv- 
ity in uniformization of the magnetic field structure dis- 
torted by the onset of the Parker instability. Referring to 
the picture of the galaxy NGC 6946 presented by Beck 
and Hoernes (jl996|l we note that highly unpolarized syn- 
chrotron emission comes from the regions of optical arms 
and a polarized emission comes from inter-arm regions. 
The reason for the depolarization in arms is obviously a 
strong agitation of the interstellar medium (ISM) by stel- 
lar explosions and winds. These phenomena form a very 
intense source of perturbations which are supposed to con- 
tribute to the excitation of the Parker and other instabili- 
ties. Since density waves are interpreted as waves traveling 
across the galactic disk area (Binney and Tremaine ,1987,) 
optical arms were placed in the past in the present lo- 
cation of inter-arm regions. One can ask therefore which 
mechanism is responsible for the converting of the tur- 



bulent medium in the arms in the much more uniform 
medium in the inter-arm regions. Within the framework 
of our model, we can suggest that the magnetic reconnec- 
tion can play an important role in the uniformization of 
galactic magnetic fields. 

2. Method 

2.1. Equations 

We investigate the Parker instability in galactic disk in 
the presence of a rigid rotation and a fast magnetic re- 
connection assuming isothermal evolution of the system 
following Paper 1. 

We solve a set of MHD equations including the Coriolis 
force in the equation of motion and the resistive term in 
the induction equation: 



dp 
at 



V-(pt;) = 0, 



(1) 



- + (. • V) . = V + - j + '-^-2n.v+g , (2) 



OB 



(3) 



where ft is the angular rotation velocity, j' = V x B is 
the current density, g = g e^, g = const is the gravita- 
tional acceleration and the other symbols have their usual 
meaning. Our assumption of uniform gravity is not realis- 
tic for actual galaxies, nevertheless it is satisfactory for a 
theoretical investigations of the Parker instability. 
We adopt an isothermal equation of state 



P 



(4) 



where Cs= const is the isothermal sound speed. 

The following Ansatz for an anomalous resistivity 
(Ugai 1992, Konz et al. 00001 Tanuma et al. .2001. Paper 
I) is applied: 



=m+ m {f - fcritf'^ 6 {f - Jcrit) > 



(5) 



where r^i and 772 are constant coefficients, jcrit denotes 
the critical current density above which the resistiv- 
ity switches on and O is the Heaviside step function. 
Such a prescription of the resistivity allows for a use- 
ful parametrization of the magnetic reconnection, as dis- 
cussed in Paper I. 

Following Paper I we introduce the Coriolis force cor- 
responding to the rotation angular velocity 17 which takes 
a typical galactic value defined in section 2.3. The Coriolis 
force is responsible for horizontal deflections of magnetic 
field lines, which are smaller in magnitude than the ver- 
tical deflections. The role of the Coriolis force is to twist 
rising magnetic loops against galactic rotation and to pro- 
duce, in cooperation with vertical gravity and resistivity, 
the helical magnetic field structures apparent in simula- 
tions of Paper I. 
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In order to integrate MHD equations JQ) - © we apply 
the ZEUS-3D code (Stone & Norman 1992a, ,1992b) with 
our own modifications necessary to introduce the Coriolis 
force and to incorporate the localized resistivity (see Paper 
I for more details). 

2.2. The method of measuring of the effect of the 
resistivity 

The presence of the resistivity in the induction equation 
releases the magnetic freezing condition and leads to a 
topological evolution of magnetic field lines. 

Measurements of deviations of a single magnetic field 
line from the perfect freezing in gas motion is possi- 
ble through integration of magnetic vector field in a 3D 
Cartesian volume bounded by periodic boundaries. We 
propose a method that is illustrated in Fig. ^ We assume 
that the initial magnetic field is parallel to the Y-axis in 
the whole volume of the computational box. We apply 
periodic boundary conditions on the XZ as well as YZ 
boundaries. We integrate magnetic field lines between the 
starting XZ boundary denoted by the symbol " S" and the 
ending XZ boundary denoted by " E" . We consider a fam- 
ily of N magnetic field lines. Each member of the family 
is denoted by the index i. Our choice of initial magnetic 
field parallel to the Y direction ensures that magnetic field 
lines are closed by periodic boundaries, i.e. the intersection 
points (Xi, Zi)s of the i-th line in the S plane are identical 
as the intersection points (Xj, Zi)E in the E plane. 

In the case of ideal MHD the magnetic freezing condi- 
tion means that a given element of the fiuid will remain 
on the same magnetic flux tube for an arbitrarily long pe- 
riod of time. This means that if (X^, Zi)s = (Xj, ZijE at 
i = then {Xi,Zi)s = {Xi,Zi)E at any later time (in- 
stant), no matter what the evolution of the box interior 
is. If, on the other hand, the system is non-ideal i.e. re- 
sistivity is present, then magnetic field lines become open 
{Xi,Zi)s ^ {Xi,Zi)E at t > 0. We define deviations of 
the intersection points of magnetic lines with the S and 
E boundaries, as (AX^, AZi) = {Xi, ZijE — {Xi, Zi)s (see 
Fig.P). 

As soon as (AX^, AZi) are known as functions of time 
one can perform a statistical analysis through the compu- 
tation of the mean and the standard dispersion of devia- 
tions of the E and S intersection points. 
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Fig. 1. Figure describing method of measuring the disper- 
sion of magnetic field lines due to resistivity. S and E are 
boundary plane at beginning and ending of y coordinate. 
Ps and Pe are points on the XZ plane where magnetic 
line starts and ends. Axi, Azi, Ari is a deviation between 
the starting and the ending points of the i-th line. The 
right panel presents the projection of the deviation on the 
XZ plane. 



In the rest of the paper, for the sake of brevity, this 
kind of dispersion will be referred to as the dispersion 
of magnetic field lines. In actual implementation of the 
method of determination of the magnetic diffusion the in- 
tegration of magnetic field lines is performed with the aid 
of the Runge-Kutta method of 4**^ order. This procedure 
is applied for a large sample of staring points Ps on the 
first XZ boundary for different timesteps of evolution of 
the Parker instability. 

The time evolution of the statistical moments depends 
primarily on values of the resistivity coefficients r\\ and 
r\i , as well as on the grid resolution adopted for the MHD 
simulations. The presented method makes it possible to 
determine the intrinsic numerical resistivity of the MHD 
code by means of comparison of the dispersion measure of 
magnetic field lines for low resolution runs with the van- 
ishing resistivity and high resolution runs with different 
values of the resistivity coefficient rji . Results of this pro- 
cedure will be presented in Section|3 Finally we would like 
to stress out that the periodic boundary conditions are es- 
sential for our method because we measure the departures 
of magnetic field lines from the closure (by periodic bound- 
aries). This measure indicates topological changes of field 
lines and therefore the presence of non-ideal effects. 



(6) 2.3. Initial conditions 



Our initial equilibrium state is an exponentially strat- 
ified isothermal disk (Parker 119661 Paper I). The ini- 
tial equilibrium is characterized by an uniform vertical 
gravity with the vertical gravitational acceleration g = 
—2 • 10~^cms^^, the midplane gas density uq — Icm"^ 
and the sound speed Cg — 7kms~^, constant across the 
disk. The fixed ratio of the magnetic pressure to the gas 
pressure is given by a = Pmag/Pgas, so that the total pres- 
sure is given by ptot — (l+ck)Pgas- We assume that a = 1.0, 
corresponding to the magnetic field strength of 4.5 juG at 
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the galactic midplane. The initial magnetic field is parallel 
to the azimuthal Y-axis. The dependence of equilibrium 
quantities on z is given by 



Table 1. Input parameters for computed models. 



no(0) "po(0) " S2(0) 



(10) 



where H = (1 + a) c^l\g\ ~ 159 pc the vertical scale- 
height. We apply the angular velocity of the galactic rota- 
tion 17 = 0.025 Myr~^ corresponding to that of the Solar 
orbit. 

In order to excite the Parker instability we apply ini- 
tial perturbations of vertical velocity component in the 
following form 




sin (^) for \z\ < H 



elsewhere 



(11) 



where rix, riy are the numbers of the harmonic compo- 
nents of the X, Y directions, respectively. Lx: Ly are the 
physical sizes of the simulation domain and the initial ve- 
locity amplitude is wq = 1.0 kms"-'^. Following Paper I the 
initial velocity perturbations are nonvanishing only up to 
the height H above the disk midplane in order to take into 
account the fact that the perturbations originate from the 
disk activity, i.e SN explosions, stellar winds, etc. For the 
sake of simplicity we perturb only the vertical velocity 
component. 

2.4. Numerical setup and input parameters 

We analyze a set of numerical simulations in a domain of 
physical sizes of — 600 pc, Ly = 1800 pc, — 1200 pc 
in the Cartesian reference frame x, y, z corresponding lo- 
cally to the radial, azimuthal and vertical galactic coor- 
dinates. Given the fact that the scaleheight of the disk is 
159 pc, the density contrast between the bottom and the 
top of the computational box is about 1900. 

The initial perturbations follow the formula with 
Tlx X riy — 3 X 3 in all cases. The periodic boundary 
conditions are applied at all boundaries perpendicular to 
the galactic plane (at the XZ planes and the YZ planes). 
At the lower and upper boundary parallel to the galactic 
plane (the XY planes) the reflecting boundary conditions 
are applied. We note that the reflecting boundary condi- 
tions imposed on horizontal boundaries are formally not 
compatible with the vertical stratification of By. However, 
this choice of boundary conditions allows for a stable evo- 
lution of the system on long timescales without any sig- 
nificant contribution of numerical artifacts. Moreover, the 
integrated magnetic field lines are located in the middle 
range of the full height of the computational box, i.e. 
the starting points of magnetic field lines span the range 
240pc < z < 960pc out of the full height of the computa- 
tional box. 

In our study we adopt several different values of the 
resistivity parameters "qi , •q2 and two different grid resolu- 
tions: 30x90x60 (low) and 60x180x120 (high) as listed in 
Tabled We assume pc^ Myr^^ as the unit of 77. 



Model 


resolution 


Vi 


V2 


Jcrit 


pO 


low 


0.0 


0.0 




Pl 


low 


0.0 


0.1 


0.1 


p2 


low 


0.0 


1.0 


0.1 


p3 


low 


0.0 


10.0 


0.1 


p4 


low 


0.0 


100.0 


0.1 


qO 


low 


0.1 


0.0 




rO 


low 


1.0 


0.0 




sO 


low 


10.0 


0.0 


_ 


to 


low 


100.0 


0.0 





PO 


high 


0.0 


0.0 




P2 


high 


0.0 


1.0 


0.1 


P3 


high 


0.0 


10.0 


0.1 


P4 


high 


0.0 


100.0 


0.1 


QO 


high 


0.1 


0.0 




RO 


high 


1.0 


0.0 




SO 


high 


10.0 


0.0 




TO 


high 


100.0 


0.0 





2.5. Statistical analysis of the simulation results 

For each model displayed in Table Q we compute a time 
evolution of the following six quantities: 

(a) the maximum of vertical velocity Wzmax representing 
the amplitude of unstable Parker modes, 

(b) the ratio of mean energy densities of total and uni- 
form magnetic fields in the disk {z < H), C/Bt/C/Bu =< 

> / < B >^. The averaging denoted by < ... > is 
performed over a volume between z — and z — H , 

(c) the maximum resistivity jymax in tbe computational 
volume, 

(d) the maximum of current density jmax in the computa- 
tional volume, 

(e) the dispersion ax of magnetic field lines in the X di- 
rection, 

(f) the dispersion cr^ of magnetic field lines in the Z direc- 
tion. 



3. Results 

All simulations presented in this paper start from the same 
initial state as described in Sect. 12.31 The basic difference 
between different considered models listed in Table fol- 
lows from the choice of different values of resistivity pa- 
rameters rji and 772- The identical initial conditions lead 
to initially identical evolution, however the differences be- 
come apparent as soon as resistivity starts to play a role. 

An example of cross-sectional slices at the end of expo- 
nential growth phase for the " ideal" model pO (771 = 772 = 
0) is shown in Fig. 13 Panel (a) displays magnetic field 
vectors on the background of density perturbation Ap/po, 
where po is the density distribution at t = 0. The initially 
horizontal magnetic field lines arc deformed in the familiar 
way by an undulatory mode of the Parker instability. Gas 
condensation of magnitude A ^ po is forming in the mag- 
netic valley. Panel (b) shows the density perturbation and 
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Fig. 2. Cross-sectional slices of the computational domain aXt = 700 Myr corresponding to the end of the exponential 
growth phase. Panel (a) shows magnetic vectors and the greyplot of density perturbation Ap/po in the YZ plane, 
panel (b) shows velocity vectors and the density perturbation in the XZ plane and panel (c) shows magnetic vectors 
and current density in the same plane. 



velocity vectors in the XZ plane. It is remarkable that the 
falling gas forms dense vertical sheets between rising mag- 
netic lobes. Panel (c) displays the XZ-components of mag- 
netic vectors superposed onto the current density greyplot. 
As it is apparent, the current sheets form between rising 
magnetic lobes and coincide with the dense vertical gas 
sheets. The cross-sectional slices for other models (except 
the cases of large values of uniform resistivity) remain very 
similar to the presented one at time T = 700 Myr. 

We refer the reader to Paper I for more a detailed dis- 
cussion of different initial configurations, which include 
cases both with and without the explicit resistivity. As 



it has been shown in Paper I, in the presence of explicit 
localized resistivity reconnection starts to operate in cur- 
rent sheets. Before the onset of reconnection, Parker insta- 
bility deforms the initially horizontal magnetic field into 
magnetic arches which are additionally twisted by the 
Coriolis force. The reconnection process transforms the 
twisted arches into helical magnetic flux tubes. The ad- 
jacent tubes undergo next phase of reconnection leading 
to their coalescence. The final states contain a significant 
radial magnetic field component, which is produced as a 
result of Parker instability, the Coriolis force and magnetic 
reconnection. 
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In the absence of explicit resistivity our system starts 
to reconnect as well, at later times, as a result of numeri- 
cal resistivity. A comparison between cross-sectional slices 
for two cases, with and without the explicit resistivity is 
presented in Figs. 2(b) and 3(a) of Paper I. 

In order to present the results of our statistical analy- 
sis of the perturbed configuration we display the six sta- 
tistical quantities (a)-(e), described in Subsection 12.51 in 
Figs. 13 - El The plots in Figs. 01 and 0] show the evolu- 
tion of these quantities for experiments listed in Table 1. 
i.e. for different values of uniform resistivity coefficient rji 
and 772 = 0. The results for low resolution runs pO, qO, 
rO, sO, to are presented in Fig. O and the results for the 
high resolution runs PO, QO, RO, SO, TO are presented in 
Fig^ Next two figures present results for different values 
of r]2 (the localized resistivity coefficient) and 771 = 0. The 
results for low resolution runs pO, pi, p2, p3, p4 are pre- 
sented in Fig. 13 and the results for high resolution runs 
PO, PI, P2, P3, P4 are presented in Fig. El 

3.1. Models with different values of the coefficient rji 

In general, all statistics are strongly dependent on the 
value of the resistivity coefficient 771 (except the maxi- 
mum resistivity, which is not representative for the present 
set of models). The evolution of Wzmax presented in 
Figs. Ela) and El^a) demonstrates a common existence of 
three phases: the initial phase, the exponential growth 
phase and the saturation one. 

During the initial phase, Wzmax decreases because the 
initial perturbation is a mixture of stable, decaying and 
unstable modes. The stable modes are quickly radiated 
out from the perturbation region, the decaying modes 
are dumped out and then the standing waves of unstable 
modes start to dominate at the end of the initial phase. 

The second phase of exponential growth, belonging to 
the linear phase of the growth of the Parker instability, 
starts at about 400 Myr and lasts about 250 Myr. This 
period is characterized by an exponential rise of the max 
by factor of about 100 (from 0.1 to 10.0 kms~^). The in- 
clination of the curve in this phase gives a good qualifica- 
tion of the evolution growth rate. For models pO, qO, rO, 
sO, PO, QO, RO, SO (771 = 0.0 ~ 10.0), the growth rate is 
about 0.036, but for the cases tO and TO with 771 — 100.0 
it is only 0.009. The maximum of v^, reaches a level of 
about 10.0 kms^^, comparable to the Alfven speed, at the 
beginning of the saturation phase and does not exceed it 
later on. For the largest value of 771 the evolution looks 
different, however we can still distinguish three phases. 
Values of i^zmax in initial phase are larger, the growth 
rate is smaller, as mentioned earlier, and in the satura- 
tion phase, the maximum velocity value decreases faster 
(but only for the low resolution) . In this case the presence 
of high resistivity apparently slows down the linear phase 
of growth of the Parker instability. 

Figs. El^b) and^fb) show the evolution of the ratio of 
mean energy densities of the total and uniform magnetic 



fields in the disk C/Bt/C^Bu for the low and high resolu- 
tion respectively. During the period of the first 500 Myr, 
the ratio is 1.0 in both cases, which means that the tur- 
bulent magnetic component is negligible. Then the ratio 
starts to rise very rapidly to the maximum value of about 
1.15 and later on it decreases. We can notice that values 
of ?7Bt/C^Bu are strongly dependent on the uniform resis- 
tivity. The experiments with the low and high number 
of the grid points show that both maximum values and 
the position of the curve maxima are different for models 
with different values of the resistivity coefficient 771 . The 
evolution of U-Bt/U-Bn in models with low resolution and 
with the lowest resistivity (pO, qO) have maximum val- 
ues of 1.18 at the position of 1100 Myr. Increasing 771 to 
1.0 (rO) and 10.0 (sO) lowers the maximum value to 1.13 
but positions of the maxima are different for these two 
cases: rO has its maximum at about 1200 Myr, while sO at 
900 Myr. We can see that there is no rule that increasing 
the resistivity causes a shift of the maximum position to 
the right or the left side, both are possible. The last ex- 
periment with the highest resistivity tO (771 = 100.0) has 
a maximum of about 1.025 more or less at the same time 
step as for the previous experiment. The high resolution 
simulations result in a slightly different evolution of the 
ratio ?7Bt/C^Bu- Due to the lower value of the numerical 
resistivity, the maximum for the model QO (771 =0.1, see 
Table ^ is slightly higher than for the case PO without 
physical resistivity at all. The maxima for both curves are 
positioned at about 950 Myr, so shifted to the left in com- 
parison with the low resolution calculations. The run RO 
shows two maxima (about 900 and 1100 Myr) and then 
the evolution is slightly different increasing its value until 
the end of the evolution. The experiments with the higher 
resistivity coefficient (SO and TO) are similar to models 
with the low number of grid points (sO, tO, Fig.El^b). This 
means that really high physical resistivity results in that 
the numerical one is insignificant. 

In case of uniform resistivity, plots of Fig. Etc) and Fig. 
^c) are presented to ensure of compatibility with the case 
of the current-dependent resistivity discussed in the next 
Subsection. 

The evolution of jmax (Fig. Eld), Efd) shows similar 
character as the ratio C/Bt/t/Bu- In the beginning the max- 
imum of the current density is very small (about 0.03), but 
after 500 Myr it starts to grow rapidly to the maximum 
value of about 0.2 and 0.5 for the low and high resolution, 
respectively. This difference is caused by the construction 
of the current density calculations, which depends on a size 
of the grid step. In time less than 200 Myr, jmax reaches 
a maximum, which is very dependent on the uniform re- 
sistivity. The evolution of the maximum value of the cur- 
rent density in computations with low grid point number 
gives a very similar value of the curve maxima (0.2) for 
all cases, except the case tO with 771 100.0 showing only 
0.03. The simulations with high resistivity show a larger 
dispersion of maximum values, three cases (PO, QO, RO) 
have maximum value between 0.4 and 0.5, while the ex- 
periment SO (771 = 10.0) has 0.27 and TO (771 = 100.0) has 
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0.03. In the saturation phase, the value of jmax decreases. 
The inclination of all curves in this phase is independent 
of the uniform resistivity for rji < 10.0 (runs pO, qO, rO), 
but for the higher one some nonlinear effects may play a 
significant role and the inclination is different. For exam- 
ple, the case sO with r]i = 10.0 (see Fig. ^d)) contains 
three sub-phases in the saturation state: from time 700 
Myr, 900 Myr and from 1250 Myr, where the inclination 
is different. 

The deviation of magnetic field lines between opposite 
XZ boundaries is presented in Figs. ISJe), IH^f) and Figs. 
0Je), 0ff) for the whole set of runs (see Table 1) corre- 
sponding to the constant resistivity. At the end of the ex- 
ponential growth of perturbations the field lines disperse 
quickly, but later on the resistivity limits further disper- 
sion. We separated the standard deviation of magnetic 
field lines into the X direction (Figs. El^e), Efe)) and 
Z direction (Figs. inif), EJf)). For both resolutions, the 
magnetic field lines are dispersed more strongly in the Z 
direction than in X direction, depending on the magni- 
tude of the coefficient rji. This can be explained by the 
dominance of the vertical gravity force over the horizon- 
tally directed Coriolis force. This relation between basic 
forces leads to vertical current sheets as displayed in Fig. 
Etc). The vertical current sheets lead to the dominating 
vertical dispersion of reconnected magnetic field lines. The 
maximum values of for experiments with zero or the 
smallest resistivity (pO, qO, PO and QO) are about 80 pc for 
the high resolution models and about 100 pc in the case of 
the low resolution runs. Values of cr^ are much larger: 250 
pc for the high resolution and about 150 pc for the low 
resolution runs. This behavior is consistent with the ex- 
pectation that lower resolution of simulations introduces 
larger numerical resistivity. The maximum values of both 
dispersions decrease strongly with increasing resistivity, 
resulting in 50 pc dispersion for the case tO and TO (resis- 
tivity 100.0, sec TableQ]). Much smaller displacements 
in the horizontal direction arc related to the fact that the 
primary force responsible for the instability is the vertical 
buoyancy force. 

3.2. Models with different values of ■i]2 

Figs. and present the time evolution of our six sta- 
tistical quantities in simulations with different values of 
r]2 (the low resolution cases: pO, pi, p2, p3, p4 and the 
high ones: PO, PI, P2, P3, P4, see Table P). Generally, 
the influence of the localized resistivity on the statistical 
quantities (a)-(f) is much weaker (except rymax) than the 
influence of the uniform resistivity. 

The time evolution of Wzmax proceeds similar to the 
case of the uniform resistivity. In the present case, how- 
ever, the evolution of the exponential growth phase is iden- 
tical for all runs. This is because the high currents, local- 
ized in current sheets, starts to exceed the critical current 
(may change for jcrit = 0.1) in the middle of the exponen- 



tial growth phase. In the saturation phase, the maximum 
of Wz falls down less than in the uniform resistivity case. 

Comparing to the uniform resistivity (Figs. Otb) and 
Ub)) the evolution of J/Bt/^u (Figs. I^b) and E^b)) 
looks different especially for the high resolution cases. This 
quantity exhibits a strong dependence on ryi but the effect 
of localized resistivity is only moderate. Now, the ratio 
C^Bt/C^Bu reaches comparable values for both grid reso- 
lutions, while in the case of the uniform resistivity the 
resolution influences the maximum values significantly. 

For the high resolution cases jmax grows to large val- 
ues, especially for smaller values of rji (Figs. Efd) and 
Etd)). For example, for rji = 0.0 and 772 = 0.0 jmax reaches 
the value of 0.2 and 0.5 for the low and high resolution, re- 
spectively. This difference we discussed in Subsection 13. II 
For a large 772 = 100.0 the evolution of jmax rises only to 
the value of 0.2 i.e. exceeding the critical current much 
less than in the case of small rj2 ■ 

The standard deviations of the magnetic field lines 
and iTz (Figs. Etc), Etf) and E^e), Eff) do not exhibit an 
apparent dependence on the resistivity at the small grid 
resolution, however at the high resolution the behavior 
is opposite with respect to the uniform resolution. Up to 
t ~ 900 Myr all three kinds of the resistivity produce the 
same dispersion of magnetic field lines. Contrary to the 
other cases we note a systematic growth of in the case 
of rj2 — 100 for the high resolution. The dispersion of mag- 
netic field lines in the X direction reaches about 100 pc 
a,t t = 1500 Myr, while ctx is only 60 pc for 772 = 0. This 
property makes a significant difference in the localized re- 
sistivity with respect to the uniform resistivity, because it 
is apparent that the localized one does not tend to elim- 
inate the initial dispersion of magnetic field lines, which 
appeared in early phases of operation of the resistivity in 
the period between t = 600 and t — 900 Myr. 

3.3. Determining of the numerical resistivity 

In our experiments we try to roughly estimate the value of 
the numerical resistivity. The Parker instability does not 
grow infinitely but is suppressed by the resistivity after 
exceeding some level, manifested by reaching the satura- 
tion phase in the evolution of our six statistical quantities 
(see Section O. The values of these quantities do not rise 
later on, but after reaching some maximum point, oscillate 
around it, or decrease exponentially. As an indicator of the 
total diffusivity in our models, we have chosen the disper- 
sion CTz of magnetic field lines defined in Section 2.2. We 
average in time the values of this quantity near maximum, 
i.e. the averaging is done for a period of about 400 Myr 
around maximum (see Figured). The computed values of 
(fz are listed together with errors Aa^ in Table El 

For the uniform resistivity models pO, qO, and rO (see 
Table El low resolution), has the same value of about 
150.0±A(7z pc. Independently of the magnitude of exphcit 
resistivity, the time evolution of is similar for the runs 
pO, qO, and rO. This fact indicates that numerical resis- 
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500 1000 1500 500 1000 1500 500 1000 1500 

time [Myr] time [Myr] time [Myr] 




Fig. 4. Evolution of statistical quantities for high resolution. Dependence on r]i. (solid - PO, dotted - QO, dashed - 
dash dot - SO, dash dot dot - TO) 
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time [Myr] time [Myr] time [Myr] 

Fig. 5. Evolution of statistical quantities for low resolution. Dependence on r}2- (solid - pO, dotted - pi, dashed 
dash dot - p3, dash dot dot - p4) 




time [Myr] time [Myr] time [Myr] 



Fig. 6. Evolution of statistical quantities for high resolution. Dependence on 772. (solid - PO, dotted - PI, dashed 
dash dot - P3, dash dot dot - P4) 
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qO: a, = 149.20 ' 20,53 
QO: a, = 1 74.56 1 29.93 




400 600 



800 1000 1200 1400 

time [Myr] 



sO: = 96.1 7 ! 24.40 
SO: f, = 1 13.28 ; 17.48 




400 600 



BOO 1000 
time [Myr] 



1200 1400 



400 500 



800 1000 
time [Myr] 



1200 1400 



Fig. 7. Averaging of standard deviation in Z direction near maximum for models with different values of rji . Thick 
lines indicate region of averaging and computed mean value of dz. Boxes around lines indicate errors, (solid - low 
resolution, dashed - high resolution) 



Table 2. Mean values of maximum of standard deviations in 
Z direction for different models averaged from 800 Myr to 1200 
Myr. 



Model 


Vi 


dz 




pO 


0.0 


142.6 


16.2 


qO 


0.1 


149.2 


20.6 


rO 


1.0 


150.0 


25.3 


sO 


10.0 


104.5 


20.4 


PO 


0.0 


192.5 


32.3 


QO 


0.1 


174.7 


30.0 


RO 


1.0 


151.4 


25.6 


SO 


10.0 


112.6 


21.0 



tivity is large compared to For i]i = 10.0 (model sO) 
the situation is different. The explicit resistivity is now 
dominant. 

Comparing all the above with the high resolution mod- 
els (PO, QO, RO, and SO) we conclude that the numerical 
diffusion in high resolution is much smaller, because dz is 
very large for the model with no explicit resistivity and 
its value decreases with growing rji. 



Although we have very few points, we can roughly de- 
termine the numerical resistivity in low resolution through 
plotting a log-hn graph of dependency dz ^ A + B \og{rii) 
for high resolution points (see Fig. O. We obtain values 
A = 146.3 and B = —31.9. Now we take the value of 
dz — 142.6 for model pO (the low resolution, no resistiv- 
ity) and we project it onto our the curve. We can read 
the value of rji = 1.3 which corresponds to dz = 142.6. 
This means that to get evolution of the Parker instability 
in the high resolution domain similar to model pO in the 
low resolution, we have to set rji = 1.3, because this value 
corresponds to numerical resistivity in low resolution. It 
is in agreement with our rough estimate above. 

We can go further and include the error Adz = 16.2 in 
our estimations. By going through the same procedure as 
described above, but now for values dz + Adz — 158.8 and 
dz — Adz — 126.4, we estimated a possible range for the 
value of the numerical resistivity, which is r\\ numerical G 
(0.4,4.2) (see Fig.lHJ. 

The procedure described above does not exclude the 
numerical resistivity in the high number of grid points ex- 
periments. After performing a set of simulations in higher 
resolutions we should be able to plot a dependency of the 
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iF, = 146.3-31 .9 log(7i,; 

'Jnomral ^ (0.4,4.2) 




Fig. 8. Interpolation of ry„umcricai- 

numerical resistivity on the resolution. Then it will be 
possible to determine the value of 771 numerical- 

4. Discussion 

The central issue of the present project is to identify the 
effects of three kinds of resistivity (uniform, current de- 
pendent and numerical) on the topological evolution of 
magnetic field lines after the onset of the Parker instabil- 
ity. We found that the numerical resistivity present in the 
currently-used code Zeus3D produces dispersion of mag- 
netic field lines, which is present in the action of a weak 
uniform resistivity. 

This finding allowed us to propose a new method of 
determination of the numerical diffusion in the numeri- 
cal simulations of the Parker instability. The method de- 
scribed in Subsect. 12.21 is based on a statistical analysis 
of the dispersion of magnetic field lines in a non-ideal 
medium. 

We measure the dispersion of a large set of magnetic 
field lines and relate this dispersion to the type and the 
magnitude of the resistivity. This kind of analysis is appli- 
cable to MHD simulations in a local approximation involv- 
ing periodic boundary conditions at least on two planes 
which are intersected simultaneously by a family of mag- 
netic field lines. 

We determine the numerical resistivity for low resolu- 
tion simulations comparing the mean dispersion of mag- 
netic field lines in the case of the vanishing physical resis- 
tivity to the case of high resolution and several values of 
the uniform resistivity. 

The evolution of our statistical quantities shows that 
for the largest values of the resistivity and higher resolu- 
tion of simulations the numerical resistivity is negligible. 
In some cases the extremely large uniform resistivity im- 
pedes the Parker instability. 



Our considerations imply that the Parker instability 
contributes to a randomization of galactic magnetic fields. 
The main sources of perturbations in real galaxies are su- 
pernovae supplying locally kinetic, thermal and cosmic ray 
energy. These perturbations are the most plausible mech- 
anism of the excitation of the Parker instability and the 
randomization of the magnetic field in spiral arms. The re- 
sistivity is then expected to contribute to the annihilation 
of the random magnetic field component in the inter-arm 
regions. In galaxies like NGC6946 and M51 the magnetic 
field appears to be much more uniform in the inter-arm 
regions then in arms (see model of Rohde et al. ll996|l . 

The modeling of the Parker instability in the present 
simulations is far from being realistic. In order to make 
our model feasible we applied simplifying assumptions 
of a plane-parallel initial magnetic field, uniform verti- 
cal gravity, isothermal evolution of the thermal g as , a 
lack of cosmic rays, a lack of selfgravity and differential 
rotation in the disk. All these circumstances imply that 
numbers and timescales of different phenomena are not 
realistic too. Among different effects we can expect that 
the presence of cosmic rays significantly shortens all the 
relevant timescales of observed phenomena, however we 
believe that future incorporation of all the complexity of 
realistic galactic disks will not change the qualitative pic- 
ture of the evolution of a magnetic field structure. 

5. Summary and conclusions 

In the present paper we considered the general problem 
of the influence of the physical and numerical resistiv- 
ity on the Parker instability. We focused on the issue of 
the dispersion of magnetic field lines (which according to 
Subsect. 2.2 is understood in this paper in the topological 
sense), resulting from three types i.e. uniform, current- 
dependent and numerical resistivities. 

We performed several experiments applying two differ- 
ent grid resolutions and different types and magnitudes of 
resistivity. We propose a new statistical method of quan- 
titative estimation of the influence of the explicit and the 
numerical resistivity on the dispersion of magnetic lines 
and the amount of a random component of the magnetic 
field in the considered MHD system. 

Our conclusions are: 

1. Except for the most extreme value of the uniform re- 
sistivity 771 = 100, the resistivity does not influence the 
growth rate of the Parker instability. For extremal val- 
ues of the resistivity, the linear growth rate of Parker 
instability is significantly diminished. 

2. The time evolution of the dispersion of magnetic field 
lines and the ratio of the total to uniform magnetic 
field is strongly dependent on the type and the mag- 
nitude of resistivity. 

3. The presence of all kinds of resistivity leads to a rapid 
dispersion of magnetic field lines at the end of the ex- 
ponential growth phase of evolution of the Parker in- 
stability. After the initial sudden dispersion of mag- 
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netic lines the subsequent evolution depends on the 
type of resistivity. 

- uniform resistivity tends to convert the disturbed 
magnetic field structure to the uniform state (uni- 
formization of the magnetic field) . This effect is clearly 
seen in the plots of the ratio of total to imiform mag- 
netic energies and in the plots of the dispersion of mag- 
netic field lines. 

- the localized resistivity diminishes the ratio of the 
total to uniform magnetic energies less efficiently than 
the uniform resistivity. On the other hand, for high 
resolution runs we observe the growth of the horizontal 
dispersion of magnetic field lines in the course of time, 
while the vertical dispersion tends to decrease. 

- the effect of the numerical resistivity is dependent on 
the grid resolution and is qualitatively similar to the 
effect of uniform resistivity. We determine the magni- 
tude of the uniform resistivity equivalent to the numer- 
ical resistivity in simulations with lower resolution. 

4. The method can be applied to compare the numerical 
magnetic diffusivity of different codes and studies of 
different physical processes e.g. turbulence. 
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